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O -1 Abstract 

5-H [ 

Effective quantum cosmology is formulated with a realistic global internal time given 
by the electric vector potential. New possibilities for the quantum behavior of space- 
time are found, and the high-density regime is shown to be very sensitive to the 
^ . specific form of state realized. 

o 
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Even in simple isotropic models, the high-curvature regime of (loop) quantum cosmology 
remains poorly understood. At high curvature one expects strong quantum effects sensitive 
to what state the universe is in, but the precise form of suitable states is unknown. (In 
this paper, we will show a new explicit example for this sensitivity.) The popular use of 



Gaussians or semiclassical states is hard to justify in this regime, and if one starts with 
a semiclassical state at low curvature and evolves toward larger curvatures, the quantum 
state depends on the dynamics in all cosmic phases passed through. Quantum ambiguities 
then prevent precise-enough knowledge of the state dynamics. 

As a second important issue, the problem of time [U EJ [3] remains unsolved, affecting 
the right choice of dynamics. The problem is usually evaded (but not solved) by using 
specific choices of global internal times which tend to be unrealistic near the big bang, 
such as a free massless scalar or dust. As part of the problem of time, it is not known how 
to transform quantum wave functions or entire Hilbert spaces between different internal 
times, and therefore results found with one choice of global internal time do not necessarily 
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hold for other choices. But if they depend on what time is used, they cannot be considered 
physical. 

As elsewhere in physics, effective equations provide a better handle on reliable predic- 
tions. For a canonical setting such as canonical quantum cosmology, such equations do not 
refer to entire wave functions but rather to moments of a state [HE]. Only a small number 
of moments, chiefly fluctuations and the covariance, is needed in semiclassical regimes, 
and as one evolves toward stronger quantum regimes, one can self-consistently check when 
higher moments become relevant. Effective techniques for quantum constraints 0, [7J [8] 
also allow the use of realistic local internal times, and one can change between different 
times by mere gauge transformations P, Dm EH]. In this article, we will not make use of 
these latter techniques because they require rather involved discussions of constrained sys- 
tems. Instead, we develop the effective framework of quantum cosmology for a new choice 
of internal time which is still global but more realistic at high density than a free massless 
scalar or dust: radiation. In this article we demonstrate that a more realistic choice of 
internal time arises from electric fieldfl 



2 Radiation Hamiltonian and effective dynamics 

The Hamiltonian constraint for spatially flat isotropic Friedmann-Lemaitre-Robertson- 
Walker models with radiation is 



with canonical gravitational variables (c, p) with \p\ = a 2 and c = — 7asgn(p) (with the 
Barbero-Immirzi parameter 7 relevant in loop quantizations [T6| ITT]) in terms of the scale 
factor a, the derivative being by proper time. We have Poisson brackets {c,p} = 8^7^/3 
while the momentum A of E, {A, E} = 1, does not appear in the Hamiltonian. The 
matter part is determined by the electric field E = \E\p assumed sufficiently small so as 
not to cause significant anisotropy. By writing the constraint in the form of a Friedmann 
equation, dividing H by |p| 3//2 , one can easily confirm that the matter term provides the 
correct behavior for radiation: the E-term amounts to an energy density p = E 2 /p 2 = 
E 2 /a 4 with E constant because H does not depend on the momentum A conjugate to 
E. (Alternatively, the Hamiltonian can be derived using the standard electric-field energy 
density: q a bE a E b /\/det q reduces to \p\E 2 /\p\ 3 ^ 2 with an isotropic spatial metric q a b = 
\p\8ab-) 

The momentum of E is the electromagnetic vector potential and would contribute a 
non-zero term to H in the presence of a magnetic field. However, a magnetic field requires 
deviations from homogeneity for the rotation of A to be non-zero. In the symmetric 



lr rhe use of electric fields in early universe cosmology has been implemented in inflationary theories 

2 The absolute value is taken with the flat Euclidean metric, \E\ = S a bE a E b , to keep the spatial metric 
q a b as a physical degree of freedom independent of E. 
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context used here, the restriction to pure electric fields is therefore meaningful. Since the 
electric field is canonically conjugate to the vector potential A, which does not appear 
in the constraint, E is constant and A can be used as a global internal time. We will 
call this choice electric time. To realize A-evolution, we follow standard techniques of 
deparameterization and solve the constraint equation for the momentum 



PA = E{c,p) = ±^— 2 \cW\p\. (2) 

As a function on the gravitational phase space (c,p), E(c,p) provides Hamiltonian equa- 
tions of motion for the classical c(A) and p(A), 

dc dp r 7-./ m 

— = {c,E{c } p)} and — = {p, E(c,p)} , 

as well as the basis for the quantum Hamiltonian of effective equations with respect to 
A. To transform equations or solutions to proper time r, we can multiply all d/d/1 by 
dA/dr = {A, H} = 2Ej^/\p\ 1 using (CQ). We confirm the correct classical equations of 
motion 



da sgn(p) dp E dp E . 8ttG Esgn(cp) c 
--77 = -{P, E\ = T\ t, Tj=r- = -sgn(p)- 



clr 2^\p~\dr p dA p V 3 ^\p~\ ' : 

substituting c for E in the last step, and 



a 



ld 2 a sgn(p) dc _ 2E dc _ 2E ^E 8ttG |c|sgn(p) 

dr 2 7a dr 7^ dA 7P ' IP V 3 v/jp[ 

87rG ^ 2 47tG . 
r = -^-(p + 3P) 
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where we have substituted E for c in the second line and used the electromagnetic expres- 
sions for energy density p and pressure P = |p to compare with the standard acceleration 
equation. (In what follows, we will set SttG/3 = 1 and 7 = 1, so that {c, p} = 1.) 



2.1 Effective dynamics 

The sign in (T5]) determines whether one considers solutions of positive or negative frequency 
with respect to time A. Without loss of generality, we will use the negative choice, such 
that E(c,p) = —\c\^/\p~\. Moreover, we can choose a definite sign of p (the orientation 
of space as measured by a triad) because we will consider only the approach to small p, 
not a possible transition from positive to negative p, or vice versa. For such a transition 
in dynamical terms to be described reliably, the Planck regime of quantum gravity would 
have to be much better understood than is possible at present. (To describe the transition 
non-singularly, we would have to refer to wave functions subject to a difference equation 
in loop quantum cosmology [TU [191 [20] -) We will work with positive p > 0. Finally, 
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since E(c,p) is a conserved quantity, the sign of c^Jp never changes dynamically and we 
can drop the absolute value, the two sign options here merging with the explicit ± in (J2]). 
Even for quantum states, the fact that the A-Hamiltonian E(c,p) and its quantization are 
conserved means that the absolute value can be dropped, provided that the expectation 
value (cy/p) is much larger than its quantum fluctuations. Truncating the whole state to a 

support of definite sign on the spectrum not just of \cy/p\ but also of c^/p then ensures that 
no opposite-sign solutions mix. (These sign issues are the same as in harmonic cosmology 
obtained with a free massless scalar [21, 22J. For more details, see these papers or [23]. For 
the construction of corresponding Hilbert spaces in deparameterized quantum cosmology, 
see [23].) 

Following the procedure of canonical effective equations [HE], the quantum Hamiltonian 
Eq is a function on the quantum phase space with coordinates given by expectation values 
and moments 

A(c a p b ) = ((c-(c)T(p-(p)) b ) symm (3) 

of a state (using totally symmetric ordering). These variables allow a Poisson structure by 
extending 

{(A)AB)} = ^r^, (4) 

defined for all operators A and B, to products of expectation values by the Leibniz rule. 
For fluctuations (Ac) 2 = A(c 2 ), (Ap) 2 = A(p 2 ) and the covariance C cp = A(cp), we have 

{(Ac) 2 ,(Ap) 2 } = 4C cp , {(Ac) 2 ,C cp } = 2(Ac) 2 , {(Ap) 2 , C cp } = -2(Ap) 2 . (5) 

Semiclassical states are defined generally by the hierarchy A(c a p 6 ) ~ 0(h^ a+b ^ 2 ) of the 
moments. This class of states is much more general than that of Gaussian wave functions, 
which would determine all moments in terms of at most two parameters. At the level of 
effective equations, sufficient generality of the states considered can thus be guaranteed, 
without giving rise to prejudices about the form of wave functions. With a semiclassical (or 
other) hierarchy, the set of infinitely many moments can be truncated to finitely many ones 
by approximation, allowing practical methods to study the approach to strong quantum 
regimes. High orders of the moments, if required, make the equations unwieldy, but the 
derivation as well as solutions of equations of motion for moments to rather high orders can 
be done with efficient computational codes [25]. Quantum corrections by higher moments 
are analogs of higher time derivatives in effective actions [26], amounting in quantum 
cosmology to important higher-curvature corrections. 

The quantum Hamiltonian is a power series in the moments of c and p, obtained by 
Taylor expanding the quantized (E) = (E[[(c) + (c — (c)), (p) + (p — (p)))) in c — (c) and 
p- (p): 

Eq : . {E} = Em , m + E i^|^A(« • (6) 

Since we assume an operator for c to exist, we will obtain the quantum Hamiltonian of a 
Wheeler-DeWitt quantization, as opposed to a loop quantization where only exponentials 
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exp(i5c) exist, but no c [27J 22] • (A modification of the classical dynamics by holonomy 
corrections would be required in the latter case.) Choosing totally symmetric ordering for 
Cy/p and expanding to quadratic terms with second-order moments, we have 



C CP l(Ap) 2 
2^p 8 p 3 



E Q = -cV-P-^ + - q ^c+..., (7) 



abbreviating c = (c) and p = (p) without risk of confusion. To this order, the Poisson 
structure provides effective equations 

dp /n HAP)2 (R) 

dc c C cp 3 (Ap) 2 

dA 2^ 4p 3 /2 16 p5/2 ^ > 

(10) 



d(Ap) 2 


(Ap) 2 


dA 


Vp 


dC cp 


l(Ap) 2 


dA 


4 p 3 / 2 


d(Ac) 2 


(Ac) 2 


dA 


VP 



(11) 

+ < 12 > 

derived from Hamiltonian equations of motion df(c,p, A(-))/dA = {/, Eg(c,p, A(-)}. 
2.2 Solutions 

The equations ([8]) and (TTUT) for p and {Ap) 2 are not coupled to the other variables and can 
be solved separately. To do so, we introduce a new evolution parameter x by dx = p~ l l 2 dA, 
so that ( fTUl) is completely decoupled: d(Ap) 2 /dx = (Ap) 2 is solved by 

(Ap) 2 (x) = {Apf e x (13) 

with initial values at x — 0. Inserting this solution in ([8]) and rewriting it for p 2 , we have 
the inhomogeneous differential equation dp 2 /dx = 2p 2 — ^(Ap)le x solved by 



j)(r) . ; 1 - l ^A(l-< «). (14) 

4 

The function is real and positive for all semiclassical values, and in fact in the whole range 
of (Ap) < 2p . If (Ap) > 2p , in which case we have to be much more careful trusting our 
effective equations but may still analyze ( !T4l) to suggest possible effects to be corroborated 
further, p(x) remains real only for x < — ln(l — 4pq/(Aj))q). 
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2.2.1 Potential effects at large fluctuations 



In this brief section we collect properties of our equations and solutions when fluctuations 
become large, keeping in mind that we would have to go to higher orders in effective equa- 
tions to justify the implications. Nevertheless, it is interesting to see what the equations 
may indicate. 

The first derivative of p by x, 



dp = eX/2 (pg-K A P)S eX ' + K A P)o 



becomes zero at 



x = In 



-pi + i(Ap)g 



which is real and finite for (Ap)g > 4pg- Turning points - - bounces or recollapses - 
therefore require strong quantum effects and large relative fluctuations. 

We can see the option for turning points directly from the equations of motion, espe- 
cially (jHJ). When fluctuations become large, the moment terms in the quantum Hamilto- 
nian provide new possibilities for potential turning points of p(A). We have dp/dA = for 
(Ap) 2 = 8p 2 , requiring large relative volume fluctuations. To test whether this point can 
be a minimum, potentially corresponding to a bounce, we compute the second derivative 
of p by A, 

d^p = l J_(A# 

dA 2 2 128 p 4 1 ' 

using the equations of motion. (All terms linear in (Ap) 2 cancel.) For (Ap) 2 = 8p 2 , 
d 2 p/dA 2 = —1/4, indicating a maximum of p and therefore a recollapse. However, if (Ap) 2 
in (|15p is significant, higher moments may easily contribute and change the behavior. The 
precise form of the turning point found here is therefore an explicit example for an effect 
that is highly sensitive to the precise form of quantum state. For other examples in terms 
of wave functions, see the solutions given in [2H1 EH] • 

Quantum fluctuations could therefore trigger a recollapse, reminiscent of the effect 
pointed out in [301 [31]. At the turning point, using our solution (111]) , we have 

Ap Q 1 
4^1-4p 2 /(Ap) 2 2 P ° 

The recollapse value p{x) may be large if (Ap)jj is close to 4pjj, and it is equal to p for 
(Ap)o = 8pq. (If (Ap)o = 8pg, x < In 2 in order for p(x) to be real.) 
We can now transform back to A by integrating 



dA = Mx = Jp^e xl2 {\\- 7^f% - e-*) dx. 

4 Po 
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The result can be expressed in terms of hypergeometric functions. For instance, using for 
illustrative purposes the value (Ap) 2 , = 8pg, we find 



A{x) = A -2 ef </2e* - 1 



+ V ( j r(|)r(|) J tHl-t)h(l-te*)\ 

(16)' 

in which V represents the Euler function. 

Unfortunately, this function and especially its inversion for x(A), to be inserted in p(x), 
are complicated. We can proceed further with additional approximations. First, as one 
example to explore the strong quantum regime, we may assume that the initial value (Ap) 2 , 
is close to 8^q. Thus x — > 0, and if we are interested in what happens close to x ~ 0, we 
can expand 

£~vk(i-^)+<V>. 

Integrating from to a; with Aq = 0, we find 

A(x) ~ ^x (l - ^x 2 ) + 0(x 4 ) , (17) 



inverted by 



At 1 / A x 2 



x=— 1 + - — +0(^ 4 



We can then go back to (JHJ), expand it as p(x) = p (l — \x 2 ) + 0(x 3 ), and find 



1 / A 2 



^^"i^H^ 1 ]+0{A 

2.2.2 Small fluctuations 

For small relative initial fluctuations (Ap)q/pI <C 1, for which our effective equations 
are reliable, and sufficiently small x, we can expand the square root and integrate dA ~ 
(! _ i((Ap)2/^)(i - e -*)) dx to 

= 2 V /^ V / T 3 6 2 sinh(x/2 + arcosh(l/ v / l - & 2 )) + A 
with 6=1 — |(Ap)o/po- To express x in terms of A, we write 

-1/2 

x = 2arsinhU(l-(l-^(Ap) 2 1 /p 2 ) 2 ) (A - A )/V^ 



2arsinh f~ (l - (1 - i(Ap) 2 ^ 2 ) 2 ) 
-arcosh ^l-(l-i(A P ) 2 /p 2 ) 2 ) " . 



We emphasize that we had to assume relative fluctuations to be small only at one time, 
which could be in a semiclassical regime. Our solutions are then valid even in stronger 
quantum regimes (until, of course, higher moments grow large). 

Finally, we can integrate all equations (l8~])-( H2p perturbatively if we assume moments 
to be small throughout the whole evolution. At zeroth order, we first solve the classical 
equations, ignoring all moments. We obtain 



with initial values po and cq when A = 0. These solutions can then be assumed in the 
equations of motion for moments to find approximate solutions for the latter. We obtain 



(Ap) 2 (A) (xp(A) , C cp (A) oc -c(A) + const and {Ac) 2 {A) oc c 4 + const'c 3 + const V . 



(The first of these equations is, of course, consistent with our full solutions f fl3|) and f fl4|) for 
small (Ap)q/j}q.) In particular, relative fluctuations (Ap) 2 /p 2 oc p -1 , C cp /(cp) oc p -1 and 
(Ac) 2 /c 2 oc const" remain small at small c. Moreover, the uncertainty product (Ac) 2 (Ap) 2 
is bounded from below by const", and the uncertainty relation will never be violated if we 
choose appropriate values for the constants. 

These solutions indicate that effective equations get better and better when one evolves 
toward larger p, but give rise to strong quantum effects with growing relative fluctuations 
at small p. 

3 Effective equations in terms of proper time 

The transformation from internal times to proper time after quantization may not be ob- 
vious because it requires a careful look at quantum corrections of different but related 
expressions — the internal-time Hamiltonian and the Hamiltonian constraint. Neverthe- 
less, once an internal time has been chosen, there is a unique procedure to transform 
equations or solutions back to proper time. We will first go through the general procedure 
to highlight ambiguities and difficulties of deparameterized quantizations, as part of the 
problem of time. 

To obtain equations of motion in proper time in our classical deparameterized model, 
we transformed d/dA to d/dr = 2E/ ^/pd/dA with constant E, computing dA/dr = 
{A, H[N]} with a Hamiltonian constraint H of lapse function N = 1. In the deparameter- 
ized quantum model, we quantize E(c,p) after having solved the Hamiltonian constraint 
equation for E; we do not quantize H itself. After deparameterized quantization, we can 
go back to a Hamiltonian constraint 




(18) 



(19) 



E 2 

VP 



E Q (c,p,A(-)) 2 _ E 2 (c^p+\C cp p 



I 2 - \{Ap) 2 cp-^ 2 + ■ • • f 



(20) 



Q 



Vp 
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with quantum corrections in Eq(c,p, A(-)). By construction, this corrected Hamiltonian 
constraint gives rise to the deparameterized model we started with if the momentum A of 
E is chosen as time. We therefore transform A-derivatives to proper-time derivatives by 
using dA/dr = {A,Hq} = 2E/^/p pa 2Eq(c,p, A(-))/y/p, the latter equation holding on 
shell when Hq = 0. 

With this procedure, the relationship between internal-time intervals dA and proper- 
time intervals dr in terms of E does not differ from the classical one, except that a new 
function Eq is used. The identification of E with Eq, without further quantum corrections, 
comes about because A has been chosen as internal time. As a time variable, it is not quan- 
tized and retains its classical form in quantum evolution equations. By construction, the 
momentum pa of A at the quantum level is then Eq, which is used in the electromagnetic 
Hamiltonian to compute the relation between A and proper time. 

If we had quantized the Hamiltonian constraint and then solved it or deparameterized it 
after quantization, additional quantum corrections with moment terms from E 2 / yfp would 
have resulted. These corrections are not included in a deparameterized quantization or 
in a quantization of the corresponding reduced phase space. A single deparameterized 
quantization is consistent as long as one does not ask how it may be related to other 
possible deparameterizations with other choices of internal time. However, if one allows 
for different internal times, the resulting quantum theories are not likely to be equivalent: 
In each case, one first solves the classical Hamiltonian constraint equation H = for 
a different variable, E in the choice made here or some other degree of freedom in a 
different model. The resulting internal-time Hamiltonians then take into account quantum 
corrections which differ from each other and from the corrections that a quantization of 
the original Hamiltonian constraint would imply. Deparameterized quantizations not only 
ignore some quantum corrections in the Hamiltonian constraint, the terms ignored even 
differ depending on which internal time one chooses. Predictions of different models are 
then unlikely to be equivalent, and a single deparameterized model cannot be considered 
physical unless one can show how its results are related to those of different time choices. 
Our considerations in this paper, as well as most results in quantum cosmology that rely on 
quantization after deparameterization, therefore cannot be considered complete. Our aim 
here is not to derive complete effects but rather to give examples for the different features 
that various choices of internal time can give rise to. 

After these cautionary remarks, we now continue with our discussion of electric time. 
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Multiplying our previous equations with 2Eq/ \/p, we obtain 



dr \ P P 2 ) \[P P 3 ^ 2 

where, consistent with our approximation, we have ignored quadratic terms in the small 
second-order moments because they would compete with higher-order moments which are 
ignored here. As one can check explicitly, the moments satisfy 

± ((A P ) 2 (Ac) 2 - C%) = (26) 

so that the uncertainty product is preserved. Any departure of a state from saturating the 
uncertainty relation 



(A P ) 2 (Ac) 2 -C 2 cp >^ (27) 



remains constant. If an initial state saturates the uncertainty relation, it will keep satu- 
rating it in this regime, and we obtain a dynamical coherent state. 

From the first of these equations it follows that the classical relation c = —\pj ^Jp = —a, 
which descends from the isotropic reduction of the Ashtekar-Barbero connection A l a = 
Y l a — 7_fQ (7 = 1) for a spatially flat cosmological model, receives quantum corrections: 
writing p in terms of the scale factor, we have 



a + \C cp /p 
c = ~ 



1 - |(A p )Vp 2 



(l + \{Ap) 2 /p^-\c cp /a 2 + ... . (28) 



(We cannot easily transform the moments to those of a and a because these variables are 
non-linearly related to p and c.) With this relation, we can use the equation for dc/dr to 
compute the acceleration equation 

ignoring products of moments. For small relative fluctuations, the equation is consistent 
with the classical one for radiation. Only the p-fluctuation enters, while all covariance 
terms cancel. Even if fluctuations could be large, they would not cause acceleration. 
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It is interesting to analyze these equations near the time when d = = p (which, as 
we recall, requires large fluctuations), where we have small c = —\C cp /a 2 and can ignore 
c 2 -terms in our equations. In addition to a, c is then nearly constant and we can integrate 
the moment equations 

d(Ap) 2 , 9 c dC cv d(Ac) 2 „c 

dr a dr dr a 

to obtain 

(Ap) 2 (r) ~ (Ap) 2 exp(-2(r-r 1 )c/a) , (Ac) 2 (r) ~ (Ac) 2 exp(2(r - n )c/a) (31) 
where t\ is the time when a = 0. 



4 Numerical solutions 

We return to the Hamiltonian system in ([7]), to seek self-consistent solutions. Here we 
only illustrate differences between classical and quantum solutions and test the validity of 
our approximations, which turn out to be very good. A detailed analysis especially near 
potential turning points could, as motivated by [32], reveal a wave function of the Universe 
that entails parity violation and non-trivial chiral effects. However, not just the numerics 
but also higher orders of our effective equations would then have to be developed. 

We attempt to solve the system (J7J) with certain initial conditions for c(A)|a=a an d the 
other variables. We may assume the normalization of p(A) initially, such that p(A)\a=a = 
1. We then fix C cp \ A =A = C° cpl (Ap) 2 A=Ao = (Ap)g, and {Ac) 2 A=Ao := (Ac) 2 ,, respecting the 
uncertainty relation (J27J). Our effective equations are valid for some time provided that 
(Ap) 2 , <C p%, and so on, but we will choose larger values just to illustrate the implications 
of quantum corrections. Even for (Ap) 2 , m p 2 ,, it turns out that we remain close to the 
classical solutions for rather long times. 

We therefore solve the system ([7]) parametrizing solutions in terms of two positive real 
numbers, namely (Ap) 2 , , (Ac) 2 , € R + , and a real number C° p G R, which together with 
p(A)\ A=Ao and c(A)\ A=Ao represent our initial conditions for the system. The moments 
must obey the uncertainty relation (l2"Tj) . With these initial values, we solve (!SI)- (fI2]) . 

Figure [1] shows good agreement with the classical solutions even for rather large fluc- 
tuations, with Ap = p and Aco = Co {C cp = 0). Numerical solutions assuming Co = 
together with (Ap) 2 , = (Ac) 2 , = 1 and = 3/4 are shown in Fig. |2j The plot for p(A) is 
increasing in A, with a certain rate amounting to the expansion of the universe. The plot 
for c(A) seems to indicate acceleration in the increasing branch, but this is realized only 
in terms of internal time A, not in terms of proper time (see (|29p ). 

In Fig. [21 we see numerical evolutions of the second-order moments. The c-variance 
decreases in A and reaches quite soon a vanishing value. On the other hand, the p- variance 
increases with the same rate of p(A). The covariance for the position and momentum 
operators is almost stable in the range studied, while the two fluctuations evolve according 
to f l3~Tj) . The good agreement with our analytical solutions and approximations is further 
illustrated in Figs. @] and [5j 



11 



1.00 



0.98 



0.96 



0.94 



0.92 



0.90 




1 ,00 



0.99 - 



0.98 - 




100 



2011 



300 



400 



500 



Figure 1: Solutions of effective equations for expectation values, plotted as their ratios to 
the classical solutions c c i assica i(v4) (top) and p c i assical (v4) (bottom). Initial fluctuations have 
been set to rather large values — Ap = Ac = Po = c = 1 with C cp = — to show 
the implications of quantum corrections more clearly. Nevertheless, the ratios stay close 
to one. 



5 Conclusions 

We have laid out the basic description of deparameterized quantum cosmology with time 
provided by the electric field. The only matter source required to formulate time evolution 
is radiation, expected to be significant in any early-universe model. No artificial matter 
sources such as dust or free massless scalars are required. 

Our discussion in this paper does not include modifications suggested by loop quantum 
cosmology, such as holonomy and inverse-triad corrections. The former would prevent us 
from using c as an operator and in moments, making the Hamiltonian more non-linear by 
replacing c with functions such as sin (5c)/ 5. Such modifications certainly alter the high- 
curvature behavior and are expected to compete with fluctuation terms. The magnitude of 
holonomy corrections depends on the parameter 8 used in this modification (a quantization 
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c(A) 




2 4 6 8 



Figure 2: The dynamics in A of the expectation values on the quantum state of the universe 
of c(A) and p(A), computed by solving the second-order equations fj8|)-f lT2|) numerically 
Here, cq = initially and C cp = 3/4, all other initial values as in Fig. [TJ 

ambiguity). However, if 5 is related to the Planck length, for instance by 5 ~ i^l \fv as 
often assumed, holonomy corrections are of the size t\jt\i with the Hubble scale i^, the 
same size expected for higher-curvature corrections. Quantum back-reaction by moments, 
on the other hand, contributes higher-derivative terms, also related to higher-curvature 
corrections. If a loop quantization is used, it is therefore impossible to study either holon- 
omy corrections or quantum back-reaction by fluctuations in isolation. In this paper, to 
provide a manageable first discussion of quantum cosmology with electric time, we have 
therefore decided to ignore loop effects so that our solutions refer to quantum back-reaction 
in Wheeler-DeWitt models. (The alternative, ignoring quantum back-reaction but keep- 
ing holonomy corrections, is not consistent unless one restricts oneself to models in which 
quantum back-reaction is absent or weak. These are only the models of harmonic cosmol- 
ogy [2TJ [33] or kinetic-dominated ones [MJ |35j [36] . ) Moreover, recent results in off-shell 
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f(A) 




1 — 1 ' 1 A 

2 4 6 8 10 

Figure 3: The variance of p (in green), of c (in blue), and the covariance between the 
two operators (in red) for a quantum state of the universe that minimizes the uncertainty 
relation (1271) . Good agreement for instance with ( 13~T1) is realized. (Initial values as in 
Fig.!) 

loop quantum gravity [37J [38], [391 HOj HTJ H2J [33] have shown that holonomy corrections in 
loop quantum cosmology must be treated with care because they trigger signature change 
at high density [H] . Their evolution equations must therefore stop before interesting high- 
density effects can be realized. 

Our discussion of the quantum dynamics was done mainly at the effective level of 
Wheeler-DeWitt quantum cosmology, allowing us to consider large classes of states with- 
out tying us down to specific wave functions (such as Gaussians). Considering quantum 
corrections by fluctuations, we have found several new possibilities for early-universe dy- 
namics, showing also the high sensitivity of potential turning points to specific forms of 
states. The small-volume behavior in electric time is rather different from that with a free 
massless scalar as time. In the latter case, the singularity is approached exponentially with 
respect to the scalar (or by a power law in proper time) while electric time, which does 
not give rise to a harmonic model, leads to a more complicated behavior. Although we 
have pointed out the possibility of inequivalent quantum corrections obtained in different 
deparameterizations (Section [3]), it remains to be seen whether these differences are solely 
due to the different matter sources or due to effects of choices of internal times. To decide 
this question, effective-constraint methods [6l [7J would be suitable but are rather compli- 
cated to perform. We therefore end this paper by concluding that much remains to be 
done before the high-curvature regime of quantum cosmology can be controlled, including 
a better understanding of the role of internal times. 
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Figure 4: The squared p- variance (Ap) 2 (A) (dashed) compared with the expectation value 
p(A) (solid). Both functions increase in nearly the same way even for large initial fluctu- 
ations as in Fig. comfirming our analytical solutions. We have (Ap) 2 > p, but relative 
fluctuations (Ap) 2 /p 2 become very small at large A. 
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